Percolation and critical O(n) loop configurations 

Chengxiang Ding 1 , Youjin Deng 2 f, Wenan GuoftL and 
Henk W.J. Blote 3 

1 Physics Department, Beijing Normal University, Beijing 100875, P. R. China 
2 Hefei National Laboratory for Physical Sciences at Microscale, Department of 
Modern Physics, University of Science and Technology of China, Hefei 230027, 
China 

3 Instituut Lorentz, Leiden University, P.O. Box 9506, 2300 RA Leiden, The 
Netherlands 

Abstract. We study a percolation problem based on critical loop configurations 
of the O(n) loop model on the honeycomb lattice. We define dual clusters as 
groups of sites on the dual triangular lattice that are not separated by a loop, 
and investigate the the bond-percolation properties of these dual clusters. The 
universal properties at the percolation threshold are argued to match those of 
Kasteleyn-Fortuin random clusters in the critical Potts model. This relation is 
checked numerically by means of cluster simulations of several O(n) models in 
the range 1 < n < 2. The simulation results include the percolation threshold for 
several values of n, as well as the universal exponents associated with bond dilution 
and the size distribution of the diluted clusters at the percolation threshold. Our 
numerical results for the exponents are in agreement with existing Coulomb gas 
results for the random-cluster model, which confirms the relation between both 
models. We discuss the renormalization flow of the bond-dilution parameter p 
as a function of n, and provide an expression that accurately describes a line 
of unstable fixed points as a function of n, corresponding with the percolation 
threshold. Furthermore, the renormalization scenario indicates the existence, in 
p versus n diagram, another line of fixed points at p = 1, which is stable with 
respect to p. 
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1. Introduction 

The reduced Hamiltonian of the 0(n) spin model is usually written as 

H/(k B T) = -J/(k B T) J2 Si ■ Si , (1) 

<ij> 

where k B is the Boltzmann constant, and T the temperature. The spin Si is an 
n dimensional vector, localized at site i of a lattice. The sum is over all nearest- 
neighbor pairs. For the special cases n = 1, 2, 3, this model corresponds to the Ising, 
the XY and the Heisenberg model, respectively. 

The 0(n) symmetry actually allows a more general form of the pair interaction- 
e.g., one may replace J/(kBT)Si ■ Sj by a(Si ■ Sj) with a an analytic function. A 
particularly useful choice [1] is a = ln(l + xSi ■ Sj), where x is a temperature-like 
variable. The so-called high-temperature graph expansion [2] for such a model yields 
a weighted sum of graphs, in which each site connects to its neighboring sites by an 
even number of 'bonds'. On lattices with coordination number of at most three, like 
the honeycomb lattice, these graphs reduce to a set of non-intersecting loops. For 
properly normalized O(n) spins, the resulting partition sum reads [1] 

Z h}op = J2x Nb n Nl , (2) 
G 

where the sum is over all possible loop configurations Q that can be constructed on 
the edges of the lattice. The number of 'bonds', i.e., edges covered by Q, is denoted 
as iV;,, and the number of loops in Q as Ni. Remarkably, the spin dimensionality n 
of the original model now appears as a continuous variable. The critical point of this 
model is exactly known for — 2 < n < 2, as well as some critical exponents [3]. 

The significance of clusters in critical or near-critical configurations of O(n) and 
g-state Potts models has already been recognized long ago [H [5j [6]. While the 
present work focuses on the percolation aspects of the two-dimensional 0(n) model, we 
consider it useful to first review the similar case of the two-dimensional Potts model. 

It is well known that the critical singularities of the ferromagnetic g-state Potts 
model can be correctly represented in terms of Kasteleyn-Fortuin (KF) clusters [5J [3] , 
also called random clusters, rather than in terms of naively defined Potts clusters. 
The latter are formed by connecting, with probability 1 , nearest- neighbor spins in the 
same Potts state, while the formation of KF clusters uses a probability p = 1 — e~ K 
instead, where K is the nearest-neighbor coupling constant of the Potts model. Thus, 
one may consider KF clusters as percolation clusters formed by a bond-percolation 
process that uses the Potts clusters as a substrate. 

Both the size of the largest Potts and that of the largest KF cluster diverge at 
the critical point, but it is the size of the KF cluster that determines the spontaneous 
magnetization, governed by the order-parameter exponent (3(q). The Potts magnetic 
correlation function appears to be equal to the probability that the correlated sites 
belong to the same KF cluster. It follows that the fractal dimension df of KF clusters 
is equal to the Potts magnetic renormalization exponent [TJ. Potts clusters are 
denser than KF clusters, and are described by a larger fractal dimension [8, 9]. At 
the Potts critical coupling K c , bond dilution of Potts clusters yields a percolation 
threshold at bond probability p = 1 — e~ Kc [10], which corresponds precisely 
with KF clusters. In the renormalization language, this KF point acts as an unstable 
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fixed point on the Potts critical line parametrized by p, and the Potts clusters are 
described by a stable fixed point at larger p. 

For the case of the tricritical g-state Potts model, which can be reached from the 
pure g-state Potts model by including vacant sites [IT] [12] , the situation is somewhat 
different. The thermodynamic singularities are still described by KF clusters, which 
are obtained by bond dilution of tricritical Potts clusters with bond probability 
p = 1 — e~ Kt , where K t is the Potts coupling at the tricritical point. But the 
percolation threshold of the bond dilution process no longer coincides with the KF 
point. The percolation threshold now occurs at p < 1 — e~ Kt , and marks an unstable 
fixed point with respect to variation of p. The fractal properties of tricritical Potts 
clusters, as well as those of KF clusters, are now described by a fixed point that is 
stable in the p direction. 

Furthermore, the pair of fixed points on the tricritical line can be related to the 
pair on the critical line. When we define q' as the number of states for which the 
tricritical Potts model has the same conformal anomaly as the q-state critical Potts 
model, it is found that the universal properties of g-state Potts clusters match those of 
g'-state tricritical KF clusters, and that g-state KF clusters correspond with diluted 
tricritical (/-state clusters at the percolation threshold [9] [13] . 

It is now natural to address a similar percolation problem defined within the 
domains separated by critical 0(n) loops. In particular, one may ask the questions for 
what bond probability p there will be a percolation threshold, what is the exponent 
associated to p, and what is the fractal dimension of the percolation clusters at the 
threshold. While it is possible to find direct answers for the honeycomb 0(1) model, 
this work will address these questions also for other values of n. In Sec. [2] we predict 
the exponent y^, which is the fractal dimension of the clusters at the percolation 
threshold, and y p , or, v p — l/y p , which controls the divergence of the percolation 
correlation length as £ = (p — p c )~ 1 / yp . Section [3] provides a numerical analysis that 
yields the percolation thresholds and a verification of the Coulomb gas result. The 
analysis is based on finite-size scaling [14] and Monte Carlo simulations of 0(n) models 
for several values of n, using a recently developed cluster algorithm [15J. In Sec. 0] we 
review the self-matching lattice argument for the case of the dual triangular lattice, 
and its consequences for our percolation problem. We conclude with a short discussion 
in Sec. 

2. Bond percolation between dual spins 

In defining the percolation problem of 0(n) loop configurations, it is convenient to 
make use of the representation of such configurations by means of Ising spins on 
the dual lattice. Spins not separated by a loop are given the same sign, and spins 
separated by one loop have opposite signs. Here we ignore the possible inconsistency 
of such an assignment with the existence of periodic boundary conditions, because 
we expect that leading critical singularities are not modified by the restriction that 
the spins are single valued. The dual Ising configuration thus completely specifies the 
loop configuration (but we do not attempt to express the 0(n) partition sum in the 
Ising language). We can now construct dual clusters by drawing, with probability 1, 
bonds between nearest-neighbor sites occupied by spins in the same state. The hulls 
of these clusters are the 0(n) loops. At criticality of the 0(n) loop model, the clusters 
are fractals with dimension d a , and the hulls are fractals with dimension di. 

On the basis of a critical loop configuration, one may now connect neighboring 
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Figure 1. An 0(n) loop configuration on the honeycomb lattice, together with 
a bond percolation configuration in dual clusters as defined within O(n) loops. 

dual Ising spins in the same state instead with a bond probability p < 1 and thus form 
new types of clusters. This is illustrated in Fig. [TJ 

The results for the percolation problem on Potts clusters described in Sec. [T] are 
relevant for the analysis of the geometric aspects of critical O(n) loop configurations, 
because of the well-known relation between critical O(n) loops and the hulls of 
tricritical random-cluster configurations [3l [T6] . As a result, the magnetic exponent of 
the tricritical Potts model is equal to the fractal dimension of the regions separated 
by loops in the corresponding 0(n) model on the surrounding lattice. Although it is 
not obvious how to interpret the 0(n) loops on the honeycomb lattice as the hulls 
of a KF random-cluster model, one may assume that the universal aspects of the 
percolation problems within the regions separated by the loops are independent of 
the details of the lattice structure. Moreover, the partition sum of the critical O(n) 
loop model on the honeycomb lattice was shown to be identical to that of a tricritical 
Potts model with vacancies |12j on the triangular lattice. This mapping identifies 
the O(n) loops with the hulls of tricritical Potts clusters. As mentioned above, the 
universal properties of the tricritical Potts clusters also apply to tricritical KF clusters. 
One can thus associate the fractal dimension all of critical O(n) loops with the hull 
fractal dimension of tricritical KF clusters, and the fractal dimension d a of the regions 
separated by loops with the fractal dimension of tricritical KF clusters. 

These exponents are exactly known by means of the Coulomb gas method 
fTF] HH [19], and verified numerically [20] for the O(n) model. The critical O(n) 
loop model and the equivalent q — n 2 -state tricritical Potts model share the Coulomb 
gas coupling g, which is given by the following equation 

<? = 4cos 2 (^), l<g<2, (3) 

Let g 1 = 1/g, thus 0.5 < g' < 1, be the Coulomb gas coupling of the g'-state critical 
Potts model, where q' is determined by 

q' = 4cos 2 ( 5 '7r). (4) 
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This (/-state critical Potts model has the same conformal anomaly as the critical 0(n) 
loop model, as well as the n 2 -state tricritical Potts model [2~ll [22] : 



c=1 _6(l-^ = 1 _6(l 7 T (5) 
9 9 

As already mentioned in Sec. [1] the Potts clusters of this g'-state critical Potts model 
are equivalent with the q-state KF clusters of the tricritical Potts model, and therefore 
also with the regions separated by critical O(n) loops. Just as the process of bond 
dilution of critical q'-statc Potts clusters leads to KF random clusters, we expect 
that bond dilution of the dual clusters defined within O(n) loops will, at the bond 
percolation threshold, lead to configurations with critical KF-like universal properties. 
The fractal dimension of the KF clusters of the critical g'-state Potts model is [TH] 

y h = l + g'/2 + 3/(8g'), (6) 

which can be continued analytically into the tricritical range [23]. This holds as well 
for the bond-dilution exponent, which is given by [24] 

y p = 1 - 3. 9 '/2 + l/(2 5 ') . (7) 

For the special values n = 0, 1 and 2 it is possible to derive exact percolation 
thresholds. For n = 0, loops are in fact forbidden due to their zero weights, 
but, depending on the boundary conditions, loop segments may emerge from the 
boundaries. In the high-temperature O(0) phase, these segments will, however, be 
confined to a boundary layer of finite thickness, and the bulk of the model will be 
empty. Therefore, the percolation threshold is exactly that of the triangular bond- 
percolation model, which is given as the solution of 

p 3 - 3p + 1 = . (8) 

Assuming continuity of the percolation threshold between the high-temperature O(n) 
phase and the critical state, the bond-percolation threshold p c (n = 0) at O(0) 
criticality is also equal to the solution of Eq. ([5]) , which is 

p c (n = 0) = 2sin(7r/18) (9) 

For the 0(1) loop model, one may apply an exact duality transformation which 
yields the critical triangular Ising model. Thus, the dual Ising configurations described 
above are precisely those of the critical triangular Ising model and we can use its known 
properties. The critical point of the triangular model [25] is K c = (1/4) In 3, and its 
random-cluster representation determines the percolation threshold as 

p c (n = 1) = 1 - exp(-2i<g = 1 - 1/V3 . (10) 

For the case n = 2 we apply an argument of a different nature. First we note that, 
along the Potts critical line as parametrized by p, the two fixed points, describing the 
Potts and KF clusters respectively, merge [9] for q — > 4, the point where the Potts 
critical and tricritical branches meet. At this point, the difference between the critical 
and tricritical KF clusters vanishes. Therefore we expect that no further bond dilution 
of the aforementioned dual clusters is required, i.e., 

p c (n = 2) = l. (11) 

The predictions for the percolation thresholds at n = 1 and 2 will be the subject of 
numerical verification in Sec. 14.31 
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3. Simulation 

3.1. Sampled variables and finite- size scaling 

The representation of the Potts model by means of KF clusters has led to the 
development of cluster Monte Carlo algorithms [26, 27], which drastically reduce the 
critical slowing down problem in simulations of the Potts model. Since then, more 
cluster algorithms have been developed, so that accurate simulation results can now 
be obtained for a considerable number of other critical model systems. Here we use an 
efficient cluster algorithm |15j for O(n) loop models with noninteger n > 1 to verify 
the predictions made in Sec. O The simulations took place at the critical point which 
is given by x c = [2 + (2 — n) 1 / 2 ] -1 / 2 for the honeycomb lattice [3j. The loop model 
cluster algorithm |15j easily allows meaningful simulations up to a linear system size 
L = 512 at the critical point. 

The configurations generated by the Monte Carlo algorithm are represented by 
means of dual Ising spins. The percolation problem involves the addition of bonds 
between equal nearest-neighbor Ising spins with probability p. For this percolation 
problem on the dual clusters, we expect, at least in part, a similar behavior as a 
function of p as usual in percolation theory 28] . Thus, for small p the percolation 
clusters are small, and they will grow with increasing p until the percolation threshold 
p c where the largest percolation cluster diverges, at least in the thermodynamic limit. 
In a finite system, the largest cluster is limited by the system size. When a percolation 
cluster reaches the system size, we call it a "spanning cluster" . However, this statement 
has to be made more precise. For a finite system with periodic boundary conditions, 
there are different rules to define a spanning cluster. One may define it as a cluster 
whose linear size in at least one of the lattice directions reaches the size of the 
periodic box, or as a cluster that connects to itself along at least one of the periodic 
directions 1291 . Here we use the latter definition. 




In order to obtain the threshold p c and some critical exponents associated 
with this percolation problem, several quantities are sampled. These include the 
susceptibility- like quantity x G i the probability R e that the occupied bonds form a 
nontrivial loop ('nontrivial' means here that the loop spans the torus and thus cannot 
be shrunk into a point by a continuous deformation), and the density of the 
spanning cluster. The quantity R e is alternatively called the wrapping or the crossing 
probability. We provide some further details to describe these quantities. 

On the triangular lattice, a cluster may span the system in different directions, 
including the ones labeled x, y and z in Fig.[2j We thus define the measure R e of the 
spanning probability as 



where (• • •) stands for ensemble averaging, and the subscript e means that a connection 
can exist along each of the x, y, and the z directions. If there is no cluster that spans 
the system, we put R x = R y = R z = 0. We put R x — 1 if there exists of a cluster that 
connects to itself over a displacement equal to a unit vector along the x direction; if 
there are no other connections to itself in other directions, we put R y = R z — 0. The 
same applies with cyclic permutations of x, y, and z. However, a cluster may also 
connect to a periodic image of itself in other directions than those of the x, y, and z 
axes. If the cluster connects to itself with a nonzero displacement vector that is not 
parallel to the x, y, or z direction, we put R x = R y = R z = 1. With this definition of 




R e = ((R x + R y + R z )/3), 



(12) 
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Figure 2. A percolation cluster may span the periodic system with hexagonal 
symmetry in the x, y and z direction. 

R e , the x, y, and z directions are treated equivalently, despite the fact that there are 
only two independent directions in two dimensions. 

Finite-size scaling yields the following finite-size behavior for R e as a function of 
the finite size L and the bond probability: 

R e = R ec + a(p - Pc )L y p +■■■ + biL* + b 2 L y 2 +■■■ , (13) 

where p c is the percolation threshold, y p is the bond-dilution exponent, and yi, y 2 , ■ ■ ■ 
are negative correction-to-scaling exponents. The quantity R ec is defined as the value 
of R e at p c , which is also universal [30l [31] , but still dependent on the finite system 
geometry and, at present, on the underlying 0(n) critical state. 

The percolation susceptibility \ G an d the percolating cluster density Poo are 
defined as 

N c 

Xg =L- d C£n*), (14) 

i 

P x = L- d { noo ) (15) 

respectively, where d — 2 represents the dimensionality of the model, and N c is the 
number of clusters. The number of sites in the i-th cluster is denoted rij, and refers 
to that in the spanning cluster. Thus Poo represents the probability that a randomly 
chosen site belongs to the spanning cluster. The definition of \ G ^ s the same as the 
Potts magnetic susceptibility, expressed in terms of an ensemble average over random 
clusters. 

Finite-size scaling predicts the following behavior for and \g '■ 

Poo = L%- d (a + ax(p ~ Pc )L y * +--- + b l L^ +b 2 L** +•••), (16) 

Xg - L 2 y h ~ d (a + ai(p - p c )L v > +■■■ + b x L^ + b 2 L^ +■••), (17) 

where yu is the fractal dimension of the percolating cluster. At the percolation 
threshold p c , these equations reduce to 

Poo = ^- d (a + bxL"! + b 2 W* +•••), (18) 
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and 

Xg = L 2 y h - d (a + hL^ + b 2 W* +•••). (19) 

3.2. Results 

We illustrate the numerical procedure, using the 0(1.5) loop model as an example. 
The model with periodic boundary conditions was simulated at its critical point [3], 
which is x c = 0.60778 • • •. Since the cluster algorithm hardly suffers from any critical 
slowing down, as described in Ref. [15] , samples were taken at intervals of only 2 
cluster steps. 

The first stage involved the determination of the spanning probability R e (p,L), 
the density P^ of the percolating cluster, and the percolation susceptibility \G for 
several values of the bond probability p. This was done for 6 system sizes L ranging 
from 8 to 256. After the equilibration of the system, 10 s samples were taken for each 
value of p in the range 8 < L < 64, and 4 x 10 7 samples in the range 64 < L < 256. 
Statistical errors were estimated by dividing each run in 1000 partial results, and 
subsequent statistical analysis. The correlations between subsequent partial results 
are negligible for the lengths of these runs. Parts of the R e data are shown in Fig. [3j 




Figure 3. Spanning probability R e versus bond probability p for various system 
sizes in the case n = 1.5. The lines connecting the data points are added only for 
illustration. All error bars are much smaller than the size of the data points. 



The number of free parameters in the finite-size scaling equations (|13p. (fTI)|) and 
(|17p makes it necessary to apply a multivariate method. We used the Levenberg- 
Marquardt least-squares algorithm, which allows nonlinear fits according to these 
equations. We thus determined the percolation threshold p c from the R e versus 
p data. The next step involved a simulation of seven system sizes in the range 
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8 < L < 512 at our estimated value of p c , with the same lengths as mentioned 
above for 8 < L < 256, and 10 7 samples for L = 512. Including these runs, we fitted 
the unknowns in the finite-size scaling formula Eq. (I13|) . also including the universal 
probability R ec and the dilution exponent y p to the data. This yielded our final 
estimates, namely p c = 0.50403 (1), R ec = 0.6420(2), and y p = 0.395 (2) for the 
bond-dilution exponent. The latter result is in a good agreement with the Coulomb 
gas prediction y p — 0.3955 ■ • ■. 

Since the analysis of the data for the density P^ of the percolating cluster 
and for the percolation susceptibility \G yielded the p dependence in terms of the 
coefficients a.;, we can deduce the value of Poo and xg at our best estimate for the 
percolation threshold. Parts of these data for P^ are shown in Fig. 2] versus L. We 
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Figure 4. Density Pea of the spanning cluster at the percolation threshold versus 
system size L for the case n = 1.5, shown on logarithmic scales. The dashed line 
represents a fit to the data points according to Eq. (1 1 8 P . The error bars are much 
smaller than the size of the data points. 

find that, at the percolation threshold, P^ and x G are wen described by power laws as 
a function of sufficiently large lattice sizes L, in agreement with the finite-size scaling 
behavior expressed by Eqs. (TT5|) and (fll?)) . A fit of the numerical data for and \g 
according to these equations yields the fractal dimension of the percolation clusters 
as yh = 1.8678 (2) and yu = 1.8678 (1) respectively. Both values are in a good 
agreement with the value 1.86775 • • • based on the Coulomb gas prediction, Eq. ([B]). 

The same procedure as described above was applied to the cases n — 1.0, 1.25, y/2, 
1.5, 1.75,1.90 and 1.95. The percolation thresholds p c , the bond-dilution percolation 
exponent y p , and the fractal dimensions of the percolation cluster yh are obtained 
similarly. The results are listed in Table [TJ The exponents are in a satisfactory 
agreement with the values predicted by Eqs. |6j) and 0. 

We also performed simulations for n = 2. However, the R e versus p data, which 
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Figure 5. Spanning probability R e versus bond probability p for various system 
sizes in the case n = 2. The lines connecting the data points are added only for 
illustration. 

are shown in Fig. El do not show intersections for different system sizes. Near p = 1, 
the curves appear to run asymptotically parallel at vanishing distances for large L. 
These findings are consistent with the predicted values p c = 1 and y p — 0. However, 
the vanishing of y p renders Eq. (|13[) insufficient for a numerical determination of p c 
from the R e data for n = 2. We thus fitted the exponent using the data for P^ 
and Xg & t the theoretical value p c = 1. These results are included in Table [1] Before 
attempting a numerical analysis of the R e data, we will adapt Eq. (|13p for the case of a 
marginally relevant bond dilution field, with the help of renormalization considerations 
in Sec. [H 

4. Self-matching argument and renormalization flow 

Jf.,1. Matching property 

First we shall briefly review the matching properties [32l [33] of planar lattices and 
the consequences for the site-percolation thresholds of some of these lattices. Let 
V = (V,A) be a planar lattice, where V is the set of lattice sites and A the set of 
edges connecting the nearest-neighbor sites of V. The faces of this lattice are polygons 
without any 'diagonals'. Let B be the set of diagonals connecting all pairs of non- 
nearest-neighboring sites within each polygon. Then, we define the lattice C, in which 
these diagonals are included, as C = (V, A + B). Then, V and C are called matching 
lattices. It is possible to define matching lattices in a more general way, but that is 
unnecessary for our present purposes. 

Suppose that a given site-percolation configuration on V percolates in the x 
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lattice direction. The existence of a percolating path prevents the existence of a 
percolating path in a conjugate percolation configuration obtained as follows. Replace 
the occupied sites of V by empty sites, and vice versa. Add the diagonals, leading to the 
lattice C, and consider the percolation problem in the other lattice direction denoted 
y. It then follows that there is no percolating path in that direction. Furthermore, 
if a given configuration of site variables on V does not lead to a percolating path in 
the x direction, then the conjugate problem must have a percolating path in the y 

(s) 

direction. As a consequence, if the site percolation threshold of lattice V is p\> , and 

(s) 

the site percolation threshold of the lattice C is p c , then the thresholds are related 
as 

pW+pW=l. (20) 

Since no 'diagonals' can be added into the triangular lattice, the triangular lattice 
is called a self-matching lattice, and the difference between the two thresholds in 
Eq. ([20]) vanishes. The matching argument thus yields that the percolation threshold 
of the triangular lattice, as well as that of other self- matching lattices, lies at p^ = 1/2. 

An important feature of the matching argument is its independence of interactions 
between the site variables, as long as these interactions are symmetric under the 
interchange of occupied and unoccupied lattice sites. 

4-2. Percolation at bond probability p = 1 

For the case n = 1 and bond probability p = 1, the model reduces to the Ising model on 
the triangular lattice with coupling K = — (In x) /2, and all neighboring sites with equal 
spins are connected by occupied bonds. The x < x c = l/\/3 region corresponds to the 
low-temperature ferromagnet, and 1 > x > x c to the high-temperature ferromagnet, 
and x > 1 to the antiferromagnet. Thus, the symmetry between + and — spins 
holds in the range x > x c . If one considers + spins as occupied sites, and — spins 
as unoccupied ones, the bond percolation model with p — 1 can be regarded as a 
correlated site percolation model with site-occupation probability p^ = 1/2. In 
particular, the x = 1, i.e., K = case reduces to the standard site percolation on 
the triangular lattice [17] . The aforementioned self-matching relation tells us that the 
whole line for p = 1 and x > critical line of the percolation type, as already 

noted in Ref. [34]. Figure 5 in the latter reference describes the renormalization flow 
of the model in the p versus K plane. Since the percolation critical line must be a 
flow line, it follows that, at the Ising critical temperature, the point p — 1 is a fixed 
point. Since the KF fixed point is unstable, and there is no sign of intermediate fixed 
points [34] , the p — 1 fixed point is stable along the p direction. 

For n / 1, the 0(n) loop model involves non-local interactions, as reflected by 
the quantity Ni in Eq. ([2]). Nevertheless, the symmetry between the + and — spins 
in the dual triangular lattice still holds as long as x > x c (n). The matching argument 
yields that, for n < 2, the point p = 1, x = x c (n) is always a fixed point for the flow 
along the p direction. 

Consider the subspace (n,p) with x = x c (n) in the three-parameter space (n,p, x). 
We have now derived two lines of fixed points as a function of n, namely p = 1 from 
the self-matching argument, and p — p c (n) from our numerical estimates. The latter 
seems to be tangent to the line n = 2, while the former is perpendicular to the n = 2 
line. This tells something about the renormalization flow near n = 2. 
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According to the aforementioned results, we conjecture the associated 
renormalization flow as shown in Fig. [6] For n — 2, the bond-dilution field is 
marginally relevant for p < 1 and marginally irrelevant p > 1. The lowest order 




Figure 6. Percolation threshold p c (n) and the conjectured renormalization flow. 
The numerical estimates of p c (n) are shown as diamonds, and the X symbols 
stand for the fixed points. The renormalization flow is represented by the arrows, 
and for n = 2 we use the thin arrows for the marginally relevant field. The curve 
is given by Eq. (1 25 I t with ao = 0.486, a\ = 1.5. 

renormalization equation in Ap = 1 — p and \J An = y/2 — n leading to the flow 
diagram sketched in Fig. [6] is 

^ = r A P VAn~ + ri (Ap) 2 (21) 

where I parametrizes the renormalization flow such that the rescaling factor is e l , and 
ro and r\ are unknown constants. The appearance of V An in this equation is in 
line with the dependence of the 0(n) critical point x c and that of the Coulomb gas 
coupling constant g on n. 



4-3. Numerical evidence for the conjectured fixed points at p{n) = 1 

If the line of stable fixed points is indeed located at p(n) = 1, x = x c (n) for < n < 2, 
the amplitude of the irrelevant bond-dilution field is zero at and only at p = 1. For 
a test, we simulated the critical n — 1.5 loop model near p = 1. The bond-dilution 
exponent is given by Eq. (JT]), where the Coulomb-gas couplings for the stable and 
unstable fixed points relate as gg 1 = 1. For n — 1.5 this yields the bond-dilution 
exponent as y' p = —0.4386 near p = 1 and y p = 0.3955 near the threshold p c . 
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Parts of the R e data near p = 1 are shown in Fig. [7] versus the bond-occupation 
probability p. The R e data lines become more and more flat when size L increases, 
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Figure 7. Spanning probability R e for n = 1.5 and near p = 1 vs. the 
bond-occupation probability p. The labels represent linear size L. The filled 
circles approximately locate at intersections of the L and the 2L data lines with 
L = 4, 8, 16, 32, 64. They rapidly approach to p = 1. 

reflecting that the bond-dilution field is irrelevant. The filled circles mark the 
intersections of the L and the 2L data lines with L = 4, 8, 16, 32 and 64. They 
are approaching p = 1 when L increases, suggesting that p = 1 is indeed a stable 
fixed point. To obtain more solid evidence, we fitted the R e data at p = 1 by 
R e (L) = R e i + bLV", and obtained R el = 0.964780 (8) and y c = -1.9 (2), which 
indicates the absence of finite-size corrections with exponent y' p = —0.4386. 

The data are also shown in Fig. [8] as \R e (L,p) — i? e i| versus L, where R e i = 
0.964780 was taken from the fit. This figure illustrates that the difference \R e (L,p) — 
R e i | at p = 1 vanishes much more rapidly than that for p =/= 1 . ft also rather clearly 
demonstrates that the finite-size correction exponent is independent of the bond- 
occupation probability p when p / 1, and that its value is in agreement with the 
expected value y' p — —0.4386. 

We observed that the R e data in range 10 < L < 512 can be well described by 



R e (p,L) = R el +a 1 (p-p cl )L v P +a 2 (p-p cl ) z L Zv p+bL- z +cL- z+v p .(22) 

The fit yields y' p = -0.40 (3), R eX = 0.96477 (3) and p cl = 0.9996 (6), consistent with 
the values found from the fit of the data at p = 1, and in good agreement with with 
the expected value p c \ = 1. 

Furthermore, we investigate whether the numerical data for n = 2 are consistent 
with the existence of a marginal fixed point at p = 1 . Integrating the renormalization 
flow, Eq. (|2"Tj) , for n = 2, and setting the finite size L equal to the rescaling factor e l , 
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Figure 8. Spanning probability, shown as \R e — R e l \ for n = 1.5 vs. system size 
L, using logarithmic scales, for several values of the bond-occupation probability 
p. Different symbols correspond with different values of p as listed in the figure. 
The value R e i = 0.964780 was taken from the least-squares fit. For comparison, 
we show a straight line with a slope y' p = —0.4386, describing the finite-size 
dependence due to a correction to scaling associated with the irrelevant field 
along the p direction. 



one finds the renormalized value Ap' as 

V = 1 A a w ■ (23) 
1-riAplnZ 

The corresponding finite-size-scaling equation for the spanning probability is 
R e (Ap,L,u) — R e (Ap' , 1, L Vi u), where u is an irrelevant field and iji its 
renormalization exponent. Substitution of Ap' , and expansion of the scaling function 
in small arguments, yields 

R e (A P ,L,u) = R ec + R x Ap' + R 2 (A P ') 2 + ■■■ + biL" (24) 

The numerical data for R e were fitted by this formula, and several variations of it, 
concerning the number of irrelevant fields, and the degree in Ap'. The data for small 
Ap reveal the existence of only one irrelevant field, with an exponent yi close to —2. 
With this exponent fixed at = — 2, we obtain fits with satisfactory residuals for the 
data with R e > 0.74. The fits deteriorate for cutoffs at smaller values, apparently 
because the expansion parameter Ap' in Eq. becomes too large. We find that 
R ec = R e i = 0.8168 (5) at the marginal fixed point for n = 2, and that it lies at 
Ap c = 1 — p c = —0.001 ± 0.002, in agreement with the expected location p c = 1. 
These data are included in Table [TJ 
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4- 4- Numerical representation of the percolation threshold 



For a description of the numerical estimates of the percolation thresholds as a function 
of n, we impose three conditions: 

(1) An oc (Ap) 2 for An -> 0; 

(2) Ap = 1/V3 for An = 1; 

(3) Ap = 1 - 2sin(7r/18) for n -> 0. 

These conditions are based on the exact values given in Sec. [2j which are supported 
and, for condition (1), supplemented by our numerical evidence. The second condition 
is a solution of /i(Ap) = 1 - 2(Ap) 2 - 3(Ap) 4 = 0, and the third one of / (Ap) 
1 - 3(Ap) 2 + (Ap) 3 =0. 

We fitted the p c (ra) data in Table I by the formula 
/(Ap) 



An = 



(25) 



A(Ap) + /(Ap) ' 

with /(Ap) = (Ap) 2 [2(Ap) 3 (4 - Ap) + f (Ap)A(Ap)] , where the amplitude A is a 
polynomial of Ap-namely, ^4(Ap) = ao + a\Ap + a-2(Ap) 2 + ■ ■ ■. This equation exactly 
reproduces /i = for An = 1 and /o = for An = 2. From our numerical data p c (n), 
we calculated the amplitude A(Ap), shown in Fig. [9j The fit for the A data yields 
ao = 0.488, ai = 1.493, and 02 = 0.000. Although ao and ai lie quite close to the 
simple fractions 1/2 and 3/2, it is very unlikely that they are equal to these values, 
according to the % 2 criterion. 
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Figure 9. Amplitude A in Eq. l|25|l . The data points are for n = 1.25, 
v2, 1.5, 1.75, 1.9, and 1.95 when Ap decreases. The straight line is drawn as 
A = a + ai Ap, with a = 0.488 and ai = 1.493. 



We conclude this section with a comment on the condition (1). For An > 0, 
Eq. (|21|) yields two fixed points, namely (Ap)o = and — V Anr^/ri. Expansion 



Percolation and critical 0(n) loop configurations 



16 



of Eq. (f2"Tj) near these two fixed points gives the bond-dilution exponents y' p = 
roV Art near (Ap) = and y p = — r ^/An, for (Ap) = — r V An/ri. Using the 
relation between g and n and Taylor-expansion of Eq. ([7]) near n = 2, one obtains 
?/p = —2^/An/n near (Ap)o = and y p = 2 V An/n near (Ap)o = — roV An/r\. 
Therefore, one has ro = — 2/7r. The amplitudes ro and ri relate to oq in Eq. (|25[) as 
ao = (ri)V(r-o) 2 . 

Table 1. Numerical results (N) for the percolation threshold p c , the bond dilution 
exponent y p , and the fractal dimension of the percolation problem on critical 
O(n) loop configurations. In the absence of intersections as in Fig. [3] the entry 
for p c at n = 2 was roughly estimated from Fig. [5] Theoretical predictions (T) are 
included where available. 



n 




Pc 


Vp 


Rec 


yjfrom P^) 


y h (from Xg ) 


1.00 


N 


0.42265 (1) 


0.542 (2) 


0.5660 (2) 


1.8749 (2) 


1.8750 (1) 




T 


0.422649 


0.5416 




15/8 


15/8 


1.25 


N 


0.45587 (2) 


0.475 (1) 


0.6018 (2) 


1.8709 (3) 


1.8710 (1) 




T 




0.4753 




1.87098 


1.87098 


V2 


N 


0.48508 (2) 


0.425 (1) 


0.6275 (2) 


1.8687 (3) 


1.8687 (1) 




T 




0.4250 




1.86875 


1.86875 


1.50 


N 


0.50403 (1) 


0.395 (2) 


0.6420 (2) 


1.8678 (2) 


1.8678 (1) 




T 




0.3955 




1.86775 


1.86775 


1.75 


N 


0.58745 (4) 


0.290 (3) 


0.6932 (2) 


1.8658 (3) 


1.8660 (1) 




T 




0.2882 




1.86603 


1.86603 


1.90 


N 


0.6924 (2) 


0.186 (5) 


0.7402 (3) 


1.8668 (4) 


1.8670 (1) 




T 




0.1882 




1.86700 


1.86700 


1.95 


N 


0.7621 (5) 


0.13 (1) 


0.7650 (5) 


1.868 (1) 


1.8684 (2) 




T 




0.1355 




1.86845 


1.86845 


2.00 


N 


1.001 (2) 


-0.01 (2) 


0.817 (1) 


1.8750 (3) 


1.87500 (1) 




T 


1 







15/8 


15/8 



5. Conclusion and Discussion 

According to the evidence presented in Sec. the dual clusters defined in the critical 
O(n) loop model should have the universal properties of KF and Potts clusters in 
the tricritical Potts model, under the condition that the O(n) loop model and the 
tricritical g-state Potts model have the same conformal anomaly, which implies that 
q = n 2 . Thus we also deduced that dilution of the dual clusters in the O(n) loop 
model by means of a bond-percolation process leads to a percolation transition with 
the same universal properties as the 'geometric' fixed point of diluted KF clusters in 
the tricritical Potts model [9j. Moreover, these universal properties should also be 
the same as KF clusters in the critical g-state Potts model, again under the condition 
that the conformal anomaly is the same. The subsequent numerical verification in 
Sec. [3] confirms these predictions in satisfactory detail. The diagram below shows the 
universal relations between the various systems by means of vertical arrows, and the 
effect of dilution is indicated by horizontal arrows. 

O(rz) dual clusters — > dilute O(n) dual clusters 

I I 
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tricritical KF clusters — > dilute tricritical clusters 

I I 

Potts critical clusters — ► critical KF clusters 
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